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Abstract 

We propose in this paper to derive and analyze a self-consistent model describing the 
diffusive transport in a nanowire. From a physical point of view, it describes the elec- 
tron transport in an ultra-scaled confined structure, taking in account the interactions of 
charged particles with phonons. The transport direction is assumed to be large compared 
to the wire section and is described by a drift-diffusion equation including effective quan- 
tities computed from a Bloch problem in the crystal lattice. The electrostatic potential 
solves a Poisson equation where the particle density couples on each energy band a two 
dimensional confinement density with the monodimensional transport density given by the 
Boltzmann statistics. On the one hand, we study the derivation of this Nanowire Drift- 
Diffusion Poisson model from a kinetic level description. On the other hand, we present 
an existence result for this model in a bounded domain. 

Keywords, drift-diffusion system, relative entropy method, diffusive limit, Hamiltonian's 
spectrum. 
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1 Introduction 

A quantum wire is an electronic component made of a periodic ion packing. The transport 
direction is large compared to the wire section, which includes only few atoms. So, the 
assumption of infinite periodic structure in the wire cross section, which allows to derive the 
usual effective mass theorem [HEO], cannot be used anymore. 

In [5], a new quantum model for nanowires is derived. Using an envelope function de- 
composition, [3] is extended to nanowires and a longitudinal effective mass model is obtained. 
However, in many applications such that FETs (Field Effect Transistors) for example, semi- 
conductor devices contains largely doped regions. In these regions, collisions play an important 
role in the transport. Usually, quantum models do not include collisions of charged particles. 
That is the reason why a diffusive model has to be developed. 
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In this paper, a Nanowire Drift-Diffusion Poisson NDDP model is derived following [7] 
by performing a diffusive limit from a sequence of Id Boltzmann equations in the transport 
direction, one for each energy band. Similarly to [S], this model takes into account the ultra- 
scaled confinement and retains information of the nanowire cross section. Moreover, a self- 
consistent model includes the resolution of the Poisson equation in the entire device. 

1.1 Nanowire quantities 

In order to define the effective masses and the other physical quantities used in the NDDP 
model, we need to recall the model derived in [5]. We consider an infinite wire defined in a 
physical domain M x u^^, where a is the typical spacing between lattice sites. The transport 
is described by a scaled Schrodinger equation in M x ojz^ containing a potential Wc generated 
by the crystal lattice, fast oscillating in the scale defined by the crystal spacing, and a slowly 
varying potential V computed self-consistently through the resolution of a Poisson equation 
in the whole domain. Since the 2d cross section lo^^ comprises few ions, Wc is considered 
periodic only in the longitudinal x-direction, also called transport direction. The variable z 
of the transverse section can be considered as fast variable, and it can be rescaled as z' = -. 
To simplify notations, we now omit the primes. Then, uj^ will denote the scaled cross section 
and we assume to work in rescaled quantities such that the periodicity is setted to 1 in the 
transport direction. 

The starting point is the definition of the generalized Bloch functions as the eigenfunctions 
of the following problem in the 3d unit cell lA = (—1/2, 1/2) x oj^ '■ 



We use here the notation y to emphasize we consider only one period (—1/2, 1/2). We point 
out that the boundary conditions are representative of our nanowire problem. Indeed, we 
consider the periodicity only in the transport direction and we choose homogeneous Dirichlet 
conditions in other directions in order to impose confinement. 

Assumption 1.1. We assume that Wc is a nonnegative potential given in L°°(JA). 

Under Assumption ll.il verified by physically relevant potentials, it is well-known that the 
eigenfunctions Xn, solutions of (jl.ip . form an orthonormal basis of Lp'iU) (see e.g. [IQ]), with 
real eigenvalues which satisfies 



Assumption 1.2. We assume that the eigenvalues En are all simple. 

Once we compute the Xn's, we can define m* which corresponds to the n*'^ band effective 
mass and which is given by 



-\^Xn + WcXn = EnXn- 

Xniy, z) = on dujz, Xn 1-periodic in y. 
lu \Xn\'^dydz = 1. 



(1.1) 



El ^ E2 ^ ... 



n— ^+00 



lim En. = +00. 



(1.2) 




(1.3) 
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where 

Pnn' = I dyXn'iy,z)xn{y,z)dydz (1.4) 



u 



are the matrix elements of the gradient operator between Bloch waves. We also define the 
quantities 

Vnn{t,x)= V{t,X,z)gnn{z)dz = {V{t,X,-),gnn), (1-5) 



with 

gnn{z) = / Xn{y,z)dy, (1.6) 
J-1/2 

where we make use of the notation (/, g) = f{z)g{z) dz. We notice that gnn's are quantities 
that contain information of the cross section and they allow to link the one dimensional 
transport direction to the entire nanowire. They play an important role not only in [5] but 
also in the NDDP model presented in this paper. 

Finally, the fully quantum longitudinal effective mass model obtained in [5] consists of a 
sequence of one dimensional device dependent Schrodinger equations, one for each band : 

tdt'lpn{t,x) = --T^dxxi^nit,x) +Vnn{t,x)'llJn{t,x). 

The physical parameters, m* and gnn-, computed once for a given device by solving the gen- 
eralized Bloch problem (jl.ip . are integrated in transport equations. Next, self-consistent 
computations include the resolution of the Poisson equation. 

Assumption 11.21 is restrictive. In [5], a more plausible assumption is discussed, related to 
symmetry properties of the crystal. Then, to each multiple eigenvalue corresponds a system of 
coupled Schrodinger equations with dimension equal to the multiplicity of the eigenvalue. The 
kinetic part of the effective mass Hamiltonian is diagonal and the coupling occurs through the 
potential. Nevertheless, the derivation of a diffusive model in this case is far from the scope 
of this paper and is not discussed here. 

1.2 Diffusive transport description 

In this paper, we mainly consider a finite wire in the transport direction x defined in the 
bounded domain oj^ = [0,L] such that L » a. We denote Q = [0, L] x this bounded 
device. Since L is large compared to a, the crystal lattice can be assumed periodic only in 
the transport direction as presented in the above subsection. Moreover, we will consider that 
the evolution of charged particles is mainly driven by collisions with phonons which represent 
lattice vibrations. A widely used model to describe such kind of transport in various area 
such as plasmas or semiconductors is the drift-diffusion equation. It consists in a conservation 
equation of the particle density in the transport direction which is called here the surface 
density Ns{t, x) and which corresponds to the integral in the direction z of the total density. 
The current is the sum of a drift term and of a diffusion term |1H [T3l [T5] . The equation reads 

dtN, - dx (^D{dxNs + Nsd.Vs)) = 0, (1.7) 
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where D is a diffusion coefficient and Vs{t,x) is the effective potential. This potential is self- 
consistant and takes into account the quantum confinement in the nanowire. Its derivation 
will be specified in Section [5] ; in particular, we will show that its expression is given by 

+00 

l^s = -log2: with 2: = ^e"^^"+^""i, (1.8) 

n=l 

where En are the eigenvalues of the problem ()1.1|) and Vnn{t,x) are the potential energies 
defined by (jl.Sp . It is also usual to introduce the Fermi level Ep{t,x) and the Slotboom 
variable u{t, x) defined by 

u = e^^ =Nse^^ = (1.9) 
Then, the current J{t, x) can be expressed as 

J = D{d^Ns + N,dc,Vs) = DNsd^Ep = De'^^'d^u. 

The electrostatic potential V(t,x,z) is solution of the Poisson equation 

-A,,,V = p. (1.10) 

The three dimensional macroscopic charge density p{t, x, z) takes into account the contribution 
of all energy bands ; it is defined as follows 

+00 

p = Y,^n9nn, (1-11) 

n=l 

where gnn{z) is given by (|1.6p and Nn{t,x) is the charge density in the transport direction 
which is expressed in the approximation of Boltzmann statistics by 

We point out that the link between the one dimensional densities A'^^ and the charge density 
p is done using the gnnS, as it is justified in [5]. We complete this system with the following 
boundary conditions 

Ns{t,x) = Nb for xeduj^, (1.13) 

V{t,x,z) = Vb{z) for xGdoj^, (1.14) 

d,V{t,x,z) = for z G dco^. (1.15) 

These boundary conditions do not correspond to the mixed type boundary conditions necessary 
for physical applications (taking in account source and drain contacts, gate...). It is chosen 
for the mathematical convenience and in particular for the elliptic regularity properties of the 
Poisson equation (|1.1U|) on our domain. 

To simplify notations, we define the functional S[V]{t,x, z) such that 

'^C^l = E ? 9nn- (1.16) 

n=l ^ 

With this notation, we have 

p = NsS[V]. (1.17) 
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1.3 Main results 

In [71[9l[19], the authors propose transport models for confined structures using the subband 
description which allows to reduce the 3d problem to a 2d transport equation. The transport 
coefficients have then to be computed by solving eigenvalue problems for the steady-state 
Schrodinger equation in the confinement direction, which is therefore one dimensional. Com- 
pared to [71[T9], the NDDP model presented in this paper, involved the resolution of the Bloch 
problem in all directions and the confinement is two-dimensional. 

The main results in this paper concern the coupled Nanowire Drift-Diffusion Poisson NDDP 
system (|1.7p - (|1.15p and are divided into two parts. In a first part we study the derivation of 
the model (|1.7p - ()1.12p in a wire with infinite extension from a kinetic level description. This 
latter model describes the interaction of charged particles with phonons at thermal equilibrium. 
Assuming that the potential V is given, we are able to state in Theorem 12.61 the convergence 
of this model towards the NDDP model. 

In a second part, we focus on the study of the NDDP model in the bounded device 
O = [0, L] X LOz- We will make the following assumptions : 

Assumption 1.3. The function D is assumed to be a function on Q and there exists two 
nonnegative constants Di and D2 such that < Di < D < D2. 

Assumption 1.4. The initial condition satisfies A'^^logA'^^ G L^iujx) and A'^^ > a.e. And 
we denote Mr = f N^dx. 

Assumption 1.5. The boundary data for the surface density Nf, is a positive constant. The 
Dirichlet boundary condition for the potential satisfies VJ, G C'^{duJxXujz) and the compatibility 
condition dzVb{z) = for all z G duJz- 

The main result is the following existence theorem : 

Theorem 1.6. LetT > 0. Under Assumptions \1.3l\1.4\ and\ 1.5l the Nanowire Drift- Diffusion 
Poisson system ()1.7p - (|1.15p admits a weak solution such that 

A.logAT, G L°°([0,r];L^K)) and yW^ G L\[0,T]; H\uj,)), 

V G L'^{[0,T];H\n)). 

To prove this result, we follow the idea proposed in [19] which relies strongly on the estimate 
on the relative entropy that we will prove in Section 13.31 The main difficulty is due to the 
quantum confinement for which we need some sharp estimates on the quantities provided by 
the Bloch problem. Such estimates are given in Section 13.11 Then a priori estimates give a 
functional framework for the quantities Ng and V. Since Ns will be defined only in LlogL, 
we need to regularized the system. For the regularized system we obtain existence of solutions 
such as in [TO] and we recover a solution of the non-regularized system by passing to the limit 
in the regularization. 

This paper is organized as follows. In Section [21 we describe the derivation of the model 
from a kinetic model taking into account the interactions of the charged particles with phonons. 
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Section[3]is devoted to the proof of Theorem ll.6l We first state estimates on the eigen-elements 
defining the Nanowire quantities. Then, we define the regularization of the model. We prove 
a priori estimates for this regularized system. Next, the regularized Nanowire Poisson system 
is analyzed. Finally, we prove Theorem 11.61 bv passing to the limit in the regularization. 



2 Diffusive limit 
2.1 Kinetic description 

The drift-diffusion model can be derived from kinetic theory when the mean free path related 
to particle interactions with a thermal bath is small compared to the system length-scale 
[U [T7]. In this section, we present the derivation of this model from the Boltzmann equation 
describing collisions of charged particles with phonons at thermal equilibrium. This equation 
governs the evolution of the distribution function fn{t,x,p) on the n^^ band whose energy is 
given by En + Vnn- Here and in the following, we shall use the notation /„ for a function 
depending on the n^^ band, and the notation / = {fn)n>i when the entire set of bands is 
taken into account. The time variable t is nonnegative, the position variable is denoted x and 
the momentum variable p. The equation writes |17| [T5] 

^t/n + ^^{vnd.rn " d^Vnndpfi) = ^QsiDn: (2.1) 

where rj is the scaled mean free path, assumed to be small. In this equation, Vn is the velocity 
given by Vn{p) = r^, the n*'^ band effective mass ()1.3p and Vnn{t, x) the effective potential 
energy associated to the n^^ band (jl.Sp . This equation is completed by the initial data denoted 
/°. 

The collision operator Qb, describing the scattering between electrons and phonons, is 
assumed in the linear BGK approximation for Boltzmann statistics. It reads 

+00 „ 

QBif)n=Yl / »n,n'{p,p)(Mnip)fn'ip)-Mn'{p)fnip))dp' (2.2) 
n'=l 

where the function 7W„ is the Maxwellian 



normalized such that 

+00 

E 

n=l 

The repartition function Z is thus given by 

+00 



Mndp = 1. (2.4) 



Z{t, x) = ^ e" (^"+^-(*'^)) . (2.5) 



n=l 



The energies En correspond to the eigenvalues of the problem (jl.ip . We notice that Assump- 
tion [TTT] allows us to give a sense to this definition of Z since Z < J2n < J2n < oo, 
where A„ are the eigenvalues of the Laplacian operator (see Section [3.ip . 
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Assumption 2.1. The cross section a is symmetric and bounded from above and below : 
3ai,a2 > 0, < ai < an^n'{p,p) < "2, yn,n' > 1, € M, \/p' G M. 

In the diffusion approximation, boundary layers appears at the frontier of the domain. 
Since the study of this phenomena is far from the scope of this paper, we consider the limit in 
the case where the spatial domain is M, assuming that there is no charged carriers at infinity 
such that lima;_5.-|-oo fnii, x, p) = 0. For the rigorous analysis of boundary layers in the diffusion 
approximation, we refer the reader to \17\ I14j. 

Let us recall an existence result for our problem. It is a direct corollary of well known 
existence results on the Boltzmann equation (see e.g. [8l[T7] and references therein). 

Theorem 2.2. Let us assume that the potential V > is given in L°°([0, T]; H'^{M. x u^)) and 
that the initial data satisfies G /^(L°°(M x M)) and > 0. For fixed i] > 0, under Assump- 
tionlKH the problem (2J^-(KW admits a unique weak solution P £ Lf^^{R+ , f {L'^ {R x M))) 
and P > 0. 

2.2 Properties of the collision operator 

We present some well known properties of the collision operator Qb defined by ()2.2p . In this 
section, the time variable t and the position variable x are considered as parameters, thus we 
omit to write the dependence on t and x. We define the weighted space 

+00 J.2 

LM = {f = (/n)n>i such that V / -{^dp < +oo}, (2.6) 

which is a Hilbert space with the scalar product 

+00 „ J. 

U^9)m=Y. / -J^dp. (2.7) 

We have the following properties for Qb whose proofs can be found in [16] (section 3.1). 

Proposition 2.3. We assume that the cross section a satisfies Assumption \2.1[ Then the 
following properties hold for Qb '■ 

(VJ:n>lkQB{f)ndp = 0. 

(a) Qb is a linear, self adjoint and negative bounded operator on L'j^. 

(Hi) Ker Qb = {f £ L^, such that 3Ns £R,fn = N^Mn} and {Ker Qb)-^ = Im Qb- 

(iv) IfPis the orthogonal projection on Ker Qb with the scalar product {■,-)m! then 

-{QB{f),f)M>ai\\f-V{f)fM. 

The third point of Proposition 12.31 implies that the equation Qb(/) = h admits a solution 
in iff /i G {Ker Qb)^- Moreover, this solution is unique if we impose / G {Ker Qb)'^ 
where {Ker Qb)^ = {/ such that Yln=i Ir fndp = 0}. As a consequence, we can define : 
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Definition 2.4. There exists O G such that for all n> 1, 



+00 „ 

QB{Q)n = -^Mn and V / Qndp = 0. (2.8) 

We define the nonnegative diffusion coefficient by 

+00 „ 

^ = E / -^'^ndp. (2.9) 
Remark 2.5. Particular case when a is constant. 

Let us assume that for all n, n', k, k' , a{n, n' , k, k') = l/r, where r is a relaxation time. After 
calculations, 

P +^ g-(i?n+Vnn) 

&n{p) = r—Mn{p) and D = t}_^ 

" n=l 2-> m=l 



We emphasize that this expression is slightly different to the formula used in [16] for example 
where D is defined as r/m. 

2.3 Asymptotic expansion for the diffusive limit 

Let us consider a solution fn of the Boltzmann equation (j2.ip and assume that it admits a 
Hilbert expansion 

fn = h,n + Vfl,n + f2,n + ••• 

Inserting this decomposition in (|2.ip and identifying with respect to powers of 77, we obtain 

QB(/o)n = 0. (2.10) 

Vndxfo,n - dxVnndpfo,n = QB{fl)n- (2.11) 
dtfo,n + Vndxfl,n " 9j,Kn<9p/l,n = QB{f2)n- (2.12) 

With ()2.10p . we get /o G Ker Qb- Thus, (iii) of Proposition [2^31 gives 

/0,n = N,Mn- (2.13) 

Injecting this expression in ()2.1ip it follows, after calculations 

QB{h)n = Hn := VnMnidxNs + N^d^Vs 



where Vg = — InZ. By Proposition 12.31 fi exists iff G {Ker Qb)^- Because we have 
/ Vn-Mndp = {vn-Mn is an odd function), this condition is true. We choose as proposed 
in Definition 12.41 Thus, 

fi,n = -en{dxN, + NsdxVs) . (2.14) 
By Proposition 12. 3| ()2.12p has a solution iff, 

+ 00 » 

/ \dth,n + Vndxfl,n " 5xKn9p/i,„ ) dp = 0. 

Using (|2.13p . (|2.14p and (|2.9p we have formally obtained the drift-diffusion equation (|1.7p . 
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2.4 A convergence proof of the derivation 

In this Section, we investigate the rigorous diffusive hmit of the Boltzmann equation (|2.ip 
as ry — 7- 0. This study is proposed in the simphfied case where the potential V is given and 
regular. The diffusive limit of a coupled Boltzmann transport equation with the Poisson 
equation is studied in [H] and with quantum confinement in [18]. The main result is the 
following theorem : 

Theorem 2.6. Let us assume that the potential V > is given in L°°{[0, T]; H'^{M. x w^)) and 
dtV is bounded in L°°{[0,T] x M x ui^) and that the initial data satisfies G /^(L°°(]R^)) and 
f > 0. Moreover, let T > and let (/^)„>i be a solution of the Boltzmann equation \2. 1]) - 
^2.5]) . Then, under Assumption \2.1[ Ng defined by Ng := Z]n>i /r Z"'^?' converges weakly 
towards Ns G L^([0,T] x M) solution of the drift- diffusion equation : 



dtNs - [D{d^N, + Nsd^Vs)) = 0, 
where Vs{t,x) = - In ( e^(^"+^""^*'^^) ) , with the initial data N^{x) = Y.n !wfn{x,p) dp. 



Notations. For the proof, we introduce the function M„ = „/ . e (am* +^"+^"") which 
such that 

VndxMn - dxVnndpMn = 0. 

Moreover, for T > 0, we consider the Banach spaces X = L°°([0, T]; and Y 



L2([0,r];L2 , 



Lemma 2.7. Assume f^ G /^(L^p) HL^^^^q^ and V is given such as in Theorem \2.6l Then, 
the unique solution f^ of the Boltzmann equations \2.1\) in L°°([0, T]; l^{L]j.^p)) is in X. More- 
over, f^ is bounded in X independently of r}. 

Proof. Assuming that all the functions are regular enough, we multiply (EH) by /^/M„ and 
integrate. We obtain 

By assumption there exists /U > such that < /.i on [0,T] x M. We define 

^'(*) = Eyy S\t) = -^Y^Jj QB{f'^)n^Jxdp. (2.16) 

Since is negative, S'^(t) > Vt G [0,r]. So, (f2T3D gives 

dX'^ 



fiX'' < -S^i. (2.17) 



dt 

Integrating this inequality allows to conclude the proof. 

To justify all calculations, we regularized the problem and consider //j a solution of the 
regularized truncated problem, f^ G V{[0,T] x M^). Thus fn satisfies ()2.15p and 

+00 „ „ ,,2 +00 „ „ ,,2 



71=1 n=l 
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Thus fn is bounded in X independently of R. We can extract a subsequence converging 
towards a function g £ X hi X-weak* . We know moreover that /r satisfies the Cauchy 
criterion in L°^([0, T], /^(L^. p)) as a solution of the truncated problem. Thus, fa converges 
strongly towards / in this space. By uniqueness of the weak* limit, g = f a.e. □ 

Lemma 2.8. Under the assumption of Lemma 2.1 , there exists f £ X and Ns G L^([0,T] xM) 

such that, up to an extraction : 

(i) P ^ f in X-weak* . 

(a) Ng Ng in Ll^-weak. 

( Hi ) If we define the current by 

+00 



J' = -T. f ''-fndp, (2.18) 



then J"^ ^ J in ^-weak. Moreover, fn = NgAdn a.e. for all n. 
Proof. We integrate ()2.17p between and t, for all t £ [0,T]. It gives 



X'^{t)-^l fx'^{s)ds+ P S\s)ds<X''{0). 
Jo Jo 



Thus, there exists a constant C > such that 



< / S'^{s)ds < C. (2.19) 







If V is the orthogonal projection on Ker Qb defined as in Proposition 12.31 there exists 
N'^{t,x) such that V{f'^)n = N'^Mn- Moreover, - Vif^) £ {Ker Qb)^ that is to say 
that Y.n>i I {fn - 'P{P)n)dp = 0. We conclude, using ([221), that 



+00 » 

N^{t,x) = Y, / f:i{t,x,p)dp:=N^{t,x). (2.20) 

We can easily show that (iv) of Proposition 12.31 is also true for the scalar product (., .)m- We 
obtain the bound 

II iP-N^Mn? ^^^^^^ < !?! r s\s) ds < Crf. (2.21) 



n=l 



Mn ai Jo 



We verify that is bounded. Thus, we can extract a subsequence satisfying (i) 

Then, from ()2.2ip . we have 



It provides that HA'^^/^Hy is bounded. Thus there exists p £Y such that, up to an extraction, 
for aU (j) £ L1l\j_^, 



/o J J ^^'^'^'^ '^'^'^^'^^ ^^Jo J J dxdpdt. (2.22) 
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For ^ G -^^^([0, T] X E) take (/>.„(t, x,p) = ^(i, x) for all n > 1, we easily verify that (p G -^^t-^M-i- 
So, if we define Ns := pZ £ L'^{[0, T] x M), we find 



This proves (ii). Moreover, from ([222]), A^I^^ = iVl'A^n ^ pMn in X-weak*. From ^Mi), 
P and NsMn have the same weak limit. So, for all n > 1, 

f = pMn = NsMn a.e. (2.23) 

Finally, we have 

= - E / ^"Z-^^' = - E / - NsMn)dp. 

Applying the Cauchy-Schwarz inequality 

With (|2.2ip . we deduce that is bounded in L},^. □ 

Proof of Theorem 12.61 Integrating the Boltzmann equation (j2.ip with respect to p, we 
find the conservation law 

dtN^, + d^J"^ = Q (2.25) 
where is the current defined in (j2.18p . Considering the function defined in ()2.8p , we have 

J'? = _V / Vnf:>dp=-{VnMn,f:i)M = —{QB{e)nj:i)M. 

The selfadjointness of Qb gives 

- J'? = i(e„, QB{nn)M = E / ®"^^y'^" dp. (2.26) 

Now, to establish the rigorous limit — )■ 0, we use the weak formulation of (|2.ip and 
(12:25]) : for all i/j £ C^{[0,T] x R^) compactly supported and for ah (j) £ C^{[0,T] x M) 
compactly supported : 

III fn{- rjdt^ - Vnd^^P + d.^Vnndp^)dxdpdt ' \^ j j j QsiDn^dxdpdt = 0, (2.27) 

and 

'''' N^dt(t>dxdt -II ,Pd^(l)dxdt = 0. (2.^ 



At this point we use Lemmas of Appendix A. On the one hand, from Lemma lA.H we can 
substitute ipi hy 4> in ()2.27p . Summing with respect to n, we immediately obtain ()2.28p . On 
the other hand. Lemmas lA.ll and IA.2I prove that we can choose -0 = 4'^^ (|2.27p for all 
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(f) G C^([0, T] X M) compactly supported. Summing with respect to n and using ()2.26p . we find 
after calculations 

' f n \ r J r, \ ^ / / / Jn 



-^E /// m[ci>-^)dxdpdt-Y, /// ^^n-^54e„0)dxdpdf 

'^^jjj "t^^i^-^nndpQn " Vn^nO:, hiZ)dxdpdt + jj ,P<l)dxdt = 0. (2.29) 

Using Lemma |2.8| we have the weak convergence of /«, Ng and J^. Moreover, we have that 
fn S L'j^ thus the limit of the first term vanishes thanks to Lemma IA.2I Moreover, this 
same Lemma proves that Vn-j^ and Vn^j^ (for s = t, x and p) are in X. Finally, since by 
assumption Vnn is bounded in L°°([0, T]; ff^(M)), we can pass to the limit and we obtain 

jj J(j)dxdt = ^ jjj ^nNsdx{@n4') - Ns(j){dxVnndpQn " Vn^ndx ln-Z)^ dxdpdt, 

which is the weak formulation of J = —D{dxNs + NsdxVg). □ 

3 Analysis of the Nanowire Drift-DifFusion-Poisson system 
3.1 Spectral properties 

In this Section, we investigate some technical Lemmas concerning spectral properties of the 
Hamiltonian defined in (jl.ip . As in Section 12.11 we denote by A„ the eigenvalues of the 
Laplacian, i.e. 

-Au„ = A„it„, on = (-1/2, 1/2) X a;^, 

n„(-l/2, z') = n„(l/2, /), My, z') = on (-1/2, 1/2) x duj,. 
Prom the min-max principle it is clear that for 

Moreover, the eigenfunctions u.y^ satisfy ,^ where (ti^)„gN are eigenvectors of the Lapla- 

cian in the y-direction with periodic boundary conditions and (n,^)„gN are eigenvectors of the 
Laplacian in the z-direction with Dirichlet boundary conditions. From well-known properties 
of eigenvalues of the Laplacian-Dirichlet [TU], we deduce that for all A > 0, J2n>i e~^^" < +oo. 
Thus 

VA>0, ^e-^-^"<oo. (3.1) 

n=l 

In the following we will make use of the notation : L^L%{Q,) = {u £ Lj^^(Q,) s.t. 111*11^^^9 (n) = 
(/^ ||u(-, z)||^,j.^ ■^dz)^/P < +oo}. We recall that we have the Sobolev embedding H^{Q) 
L'lLf{Vt) where = cj^ x C M x (gee [6]). 

Lemma 3.1. Under Assumption \l.l\ we have that for all n £ N 

IIXnllLoo(f^) < CiEn + \\Wc\\l^), (3.2) 

where C stands for a nonnegative constant. Therefore we have 

hnnh^i^^) < C{En + || || L- )' , (3.3) 
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\Vr 



nn\\L°° 



{..)<C\\V\\H^n){En + \\Wc\\L^f. 



(3.4) 



Proof. We notice first that since Wc (and Xn) is 1-periodic in x, we have that || W£||j;^oo(^) 
ll^£||L°°(n)- From we deduce that Axn G iv^(il) and 



||Ax„||i2 <2(^„ + ||Ty^||L^). 



Using the elhptic regularity for the Laplacian operator, we deduce ()3.2p thanks to the Sobolev 
embedding H'^{Q,) L°°(r2). The estimates (|3.3|) and (|3.4|) follows directly from a Cauchy- 
Schwarz inequality and the Sobolev embedding H^{Q,) ^ L^L^(J7). □ 

Lemma 3.2. Properties of S[V]. Let us assume that V € H^{Q) with V > 0. Then the 
function S[V], defined in il.l6\} . satisfies : 
(i) S\y\ is hounded i.e. \S\V]\ < +oo. 

(a) If moreover V £ H^{Q,) with V > 0, then there exists a nonnegative constant C such that : 



\\S[V]-S[V]\\L2(n)<C\\V-V\\L2^^y 
Proof, (i) First, we study the coefficient Z. We have 



(3.5) 



n=l 

Then from Lemma l3.H when V G H^{Q), there exists a constant C > such that Z{x) > C. 
Using the fact that Vnn ^ when y > 0, we get 

+ 00 +00 

|5[y](x,z)| < C^e-^'i<7„„,||L-(a..) < Cj^e-^"(£;„ + ||VFz:||l-)2, 

n=l n=l 

where we use (j3.3p for the last inequality. A direct consequence of (|3.ip is that X]„>i E^e~^" 
is finite. Then S[V] is bounded, 
(ii) We use the fact that 

S[V]-S[V]= [ dsS[V + s{V -V)]ds. 
Jo 

We define £(3) = e-(^"+<^+^(^-^)'S-)) . So, 

En>l'^'(«)fnn) ( E„>1'^(«)) " ( En>l'^(s)5r 



dsS[v + siv - V)] 

We have £'{s) = —{V — V ,gnn)£{s). The first term becomes 



En>i^'(^; 



Finally, 



En>lf(^) 



En>l'^'(^)5n 



+00 



< CY,e-''-gnn\{V -V,gnn)V 



n=l 



En>l^(«) 



+00 

< C\\V - V\\L2^^^^^)^e~^''gnn\\gnn\\L^{uj,) 
n=l 



< C\\V -v\ 
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where we use ()3.3p and ()3.ip for the last inequality. We can treat the second term in a similar 
way 



En>l ^ is)9nn En>l £{s) <V -V, gnn > 



Consequently, we deduce (|3.5p . 



<C^e-^"|(y-y,g™)|. 



n=l 



□ 



3.2 Regularized system 

We define the linear regularization operator by 



(3.6) 



where V is the extension of V by zero outside f2 and and (,e,z are C°° nonnegative compactly 
supported even approximations of the unity, respectively on M and M^. We can prove the 
following properties, using convolution results : 

Lemma 3.3. Properties of TZ"^ : 

(i) TZ'' is a bounded operator on LxL'i{Q) for 1 < p,q < +oo and satisfies for all V G L^L%{Vl), 
W[n\LlL%m < my^LUn) and Inn \\n'[V] - VhPLtin) = 0- 

(a) TZ'' is self adjoint on L^(il). 
(Hi) For allV € H^n), 

V^TZ'iV] = 7^'[V^V] and IrniWV^TZ'iV] - V^VWi^i^n) = 0- 

Then the regularized Nanowire Drift-Diffusion Poisson NDDP^ system is defined for e G [0, 1] 
by 

9iiV| - (D{d,NI + Nld.V:)) = (3.7) 



and 



TZ' 



n=l 



where the regularized quantities are defined by 



/ TZ'[V'ix,z)]gnn{z)dz =< TZ'[V%gnn >, 

+00 

log with Z' [V] = ^e- 1^"+^"" J , 



n=l 



and 



(3.8) 

(3.9) 
(3.10) 

(3.11) 



n>l 
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As above, we denote = u''e~^^"-~^^™^^ where u*^ is the Slotboom variable u"^ = Ng/Z''. 
The initial regularized density Np^ is chosen such that Ng'^ = min(A'^^, e""*^). Moreover, the 
regularized boundary conditions are 

NI{t,x) = Nb for xGduj^, (3.12) 

V%t,x,z) = Vbiz) for xGduj^, (3.13) 

d^V'{t,x,z) = for z G du^. (3.14) 

Remark 3.4. When e — 0, 7^*^ — Id and the regularized system (|1.7p - (|1.15|) tends to the 
unregularized problem p.7p - ()3.14p . 

3.3 A priori estimates 

Let us consider a weak solution (iV|, V^) of the regularized problem (|3.7p - ()3.1ip . We introduce 
two extensions Ng and y_ of the boundary data. These extensions are respectively defined on 
Ux and and chosen such that : 

• Ns_ G C'^{ujx), < Ni_ < Ns{x) < N2 where Ni_ and Ng are two constants, and Ns\q^^ = 
Nb. 

• V_£ C^(r2) and satisfies the boundary conditions Vjd^^xLj^ = ^b(-2) and ^zV^lw^x^w^ = 0. 

For regular enough domains, these functions exist. From (|1.5p with y_ instead of V", we define 
T4n[ll] denoted by Vnn - In the same way, we denote Z_, S, Nn, p, u and Ep the quantities 
associated to Ng and V_. 

Proposition 3.5. Let T > and e G [0, 1]. Let (N^, V'') be a weak solution of the regularized 
system NDDP, ^-^3A^, such thatNllogNI G L^{[0,T], L^u^)), V G L°^([0, T], //^(f])) 
and ^ Ng G L'^{[0,T], (ux)) ■ Then, there exists a nonnegative constant Ct depending only 
on T such that 

VtG[0,r], ^<W{t) <Ct, (3.15) 

where W is the relative entropy defined by 

^ = E/ {K.^n{^)-K + I^)dx + U \V{V^-V)\'dxdz. (3.16) 

Proof. We remark that 

if^i (f ) - K+N^y- = E / B,N'„i. {^y.. 

n=l ^ — n=l ^ — 

By definition, we have InA''^ = Inu'^ — En — Vnn- Using (j3.7p . it leads to 

+00 



Y,J aiKln(^)dx = j dx{De-''''dxu')\Ti(^)dx 

^E/ K{ynn-yJ)dx + Y, ^^^(F^-Kn)^^. 



n=l-"^^ n=l' 
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Integrating by parts, the first right hand side term gives 



dx. 



u 



Using the definition ()3.9p of the last right hand side term gives 

+00 „ +00 „ 

/ KdtiVL - Vnn)dx = Y, Kdt < n'[v'] - n^mann > dx. 

n=l "^"^ n=l "^"^^ 

At this point, the hnearity and the selfadjointness of the regularization operator TZ'' and the 
regularized Poisson equation ()3.8p imply 

V / KdtTZ'iV' -V]gnndxdz = -- / \V {V' - V)\^ dxdz . 

In the same way, we can write 

1 E / ^niVnn " Kn)^^ = jj^ \^ i^' " V)\'dxdz. 

Thus, defining W as in ()3.16p . we finally find 

^ = -f De-^^'^-^dx + / De~^^'^^^dx. (3.17) 



V^t) = I De~^''^^^^dx (3.18) 



We denote 

the term which can be seen as an entropy dissipation rate. We also define (3 = \\dxu/u\\L°°{u]^), 
j3 < +00. Consequently, 



Using the inequality 2ah < k^o^ + ^ for k > small enough and Assumption 11.31 we get, 

^5^<^!5^ + C,I..<C.||«,.||,,„,,. (3.19) 

where Ci and C2 are two nonnegative constants. Finally, using the inequality ln(x) — a: + 1 > 
X + (1 — e), for X > 0, we have 

Nn(j;^ + l-e)dx> N'^dx - {e - I) Nsdx. 

n=l iZJl Ji^x •'OJx 

With (|3T9]) . it leads to 

^<C2\\N!Ui^^^)<C{W + Co), 

where C and Co are two nonnegative constants. We conclude thanks to a Gronwall's inequality 
and the fact that Assumption 11.41 and V G H^{Q) imply that the initial entropy VF(0) is 
bounded. Moreover, we get the bound on the mass 

VtG[0,r], / Nldx<C. (3.20) 

□ 
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Corollary 3.6. Let T > and e G [0, 1]. Under assumptions of Proposition \3.5\ there exist 
Ci and C2 two nonnegative constants such that 

VtG[0,T], / / \d^^/Nl\^dxds < Ci, (3.21) 

Jo J U)x 

yt G [o,r], yp g [i,+oo) /* WNiUp^^^^ds < C2. (3.22) 

Jo 

Proof. In this proof, the letter C is used to denote nonnegative constants. We can express the 
coefficient defined in ()3.18p in terms of Nl and 

v'{t) = [ oUd^y^ii^ + 2d^Nid,v: + Ni\d,v,f)dx. 

J UJ.J. 

After an integration by parts on the second term of the right hand side, we deduce 

4||5xV^|li2(^^) < c[v\t) + 2^ Kd^^V^dx - 2(Nl{L)d,V^{L) + iV|(0)9,y/(0))^ . 

(3.23) 

On the one hand, after calculations, we find 

^^^^s - J z^ ' 

By the Cauchy-Schwarz inequality, the sum of the last two terms is nonpositive. Moreover, 
from the regularized Poisson equation (j3.8p 



9xxKn =< dxxTl' [V] ,gnn> = < " A,7^^ [V] , Qnn > - < W [p'] , > 

< ||'7^^[^^]||L2(^^)||Xn||L°°(t7)IIXn||//2(n)- 

To obtain the last inequality, we remark that the second term is nonpositive and to treat the 
first term, we make an integration by parts and we use the fact that 

\\^zgnn\\L'^{u)^) < 2 ||Xn (H) || Xn || ^2 (C) • 

Using the property (i) of Lemma [3^31 and (|3.ip . we conclude that 

/ Nld^^V^dx<C\\Nl\\:^.^^^)\\V^\\L2^^ay (3.24) 

J IJJx 

On the other hand, we have 

Z_jn 6 J u)z^ 

Thanks to the trace Theorem and Lemma |3.3| we obtain 

where we use the elliptic regularity and Lemma [32] and [33] for the last inequality. With (|3.23p 
and (j3.24p . we conclude that 
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Applying the Gagliardo-Nirenberg inequality to the function using the bound on 

IIA^IIIlik) ^M, we obtain 

VtG[0,r], / Id^yHVll^dx <C{1+V'{t)). 

With (fXT9]) . we can say that ^^V^{s)ds < C{W{<S) +Nit) for all t £ [0,T] and consequently 
we obtain (|3.21|) . Finally, (|3.22|) is a consequence of (|3.2ip with the Gagliardo-Nirenberg 
inequality. □ 

3.4 Analysis of the regularized Nanowire Poisson system 

In this section, the surface density Ng is assumed to be given and we only consider the 
resolution of the regularized Nanowire Poisson equation p.Sp with boundary conditions ()3.13p - 
for e G [0,1]. 
We introduce the functional space 

= {V e H^{n), s. t. Vx G du^, Vz G uj„ V{x,z) = 0}. 

Let us also take Vq G H^{Q) such that Vq = Vb on diOx x Uz and dzVo{x, z) = for all z G dooz- 
A possibility is to take Vo = y_. Most of the results presented here can be obtained by a 
straightforward adaptation of [3 [19]. Thus we will not detail the proofs. We first state the 
following existence result : 

Proposition 3.7. Let T > and e G [0,1]. We assume G L°°([0, T]; ^^(w^)) such that 
Ng > a.e. Then the regularized Nanowire Poisson equation ()3.8p with boundary conditions 
p.l3p - ()3.14p admits a unique solution G Vq + H^^ with a bound independent of e. 

Proof. Using the selfadjointness of the regularization operator, a weak solution of ()3.8p is a 
critical point in the space Vq + H^^ of the functional 

J{V,Ns) = MV) + Ji{V,Ns) 

= -f\VV\^dxdz+f iV,ln(^e"(^"+<^^[^l*->))dx. (3.25) 



n=l 



Following the proof of Proposition 3.1 in [19] and using Lemma 13.11 we show that J is a 
continuous, convex and coercive functional on Vq + H^^- Thus J admits a unique minimizer 
V and we have a bound on V in only depending on the norm of Ng. □ 

Then, we have the following continuity result : 

Proposition 3.8. LetT > and e G [0,1]. Assume Ng and Ng are given in L°°{[0,T]; L^{lOx)) 
such that Ng > and Ng > a.e. Then, the corresponding solutions and of the 
regularized Nanowire Poisson equation ()3.8p with boundary conditions ()3.13p - ()3.14p verify 

Vt G [0,r], \\V^ - V'Wmin) < C\\Ng - Nghi^^^y (3.26) 
Moreover, if Ng and Ng belongs to L'^{[0,T]; L'^{uJx)), we have 

yt G [0,r], \\V^ - V^Wn^^n) < C\\Ng - iV,||i2(,^), (3.27) 
where C stands for a nonnegative constant not depending on e. 
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Proof. Multiplying the regularized Poisson equation ()3.8p by — and integrating, we 
obtain 

j \y{V' -V')\'^dxdz = j {^{Ns-Ns)S^[V'] + Ns{S'[V']-S^[V'])^TV[V' -V']dxdz. 

Because the functional F i— )• is decreasing with respect to V , term 5^[y^])7^'^[l/'^ — 

y^] is nonpositive. We deduce 

/ \V{y' -V')\^dxdz < \\Ns-Ns\\Li^^J\<n'[V' -V%S'[V']>\\Loo^^^y (3.28) 
Jn 

Then, (i) of Lemma 13.21 and a Cauchy-Schwarz inequality gives 

|v(F^ - v^)\^dxdz < c\\Ns - Ns\\mu.jmv' - V^LUTm- 

We use the property (i) of Lemma [331 and the embedding ^ L^L^ . We obtain 

I |v(y^ - V^)\^dxdz < C\\Ns - Nshi^^jV' - V'Wmin)- 
Jn 

Finally, thanks to the Poincare inequality, we get (j3.26p . 
For the estimate, we have 

Then we bound the norm of the right hand side as above using the spectral properties in 
Section 13.11 We finally get the estimate p.27p from the elliptic regularity. 

□ 



Finally, a straightforward adaptation of Proposition 3.2 of [19] gives the following conver- 
gence result as e goes to 0. 

Proposition 3.9. As e ^ 0, the solution of the regularized Nanowire Poisson system 
converges, uniformly with respect to Ng G L^{[0,T]; L^{ojx)) such that Ns > a.e, to the 
solution V of the unregularized problem in L°°{\fd,T]\H^{^)). 

3.5 Existence of solutions for the regularized system 

Proposition 3.10. Let T > and e G [0, 1] he fixed. Then the regularized problem NDDP^ 
admits a unique solution {NI,V^) with £ C{[0,T]; L'^iuj^)) n L^{[0,T]; {u;^)) and G 
L°°{[0,T];H^n)). 

Proof. The proof of this result follows closely the proof of Theorem 1.2 of [7]. Thus we will 
not detail it and only give the main steps. The proof relies on a fixed point argument on the 
map F defined by : 

Step 1 : For a given > 0, we solve the regularized Nanowire Poisson equation (j3.8p with 
boundary conditions (f3T3]) - ([33il) and we define V"/ by ^7iO\i which belongs to L°°([0, T]; (uj^)). 
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Step 2 : The effective potential Vg being known, we solve the following drift-diffusion equation 
for the unknown 

dtNl - (^D{d^NI + iV|9,F/)) = 0, 

with the boundary condition N^Iq^^ = Nh and the initial value iV|(0, x) = Ng{x). The map 
F is then defined after these two steps by F{Ng) = iV|. 

Then we can prove that for T small enough, F is a contraction on the space Ma^T defined 
by Ma^T = {n, \\n\\T < a} where the norm is 

rT 



max^||n(t)||2,(^^) + ^ \\n{t)\\l,^^^^dtf^\ (3.29) 



.0<t<T 

We have then constructed a unique solution on a small time interval [0, Tq]. Using the a priori 
estimate, we can iterate this procedure to construct a solution on [To,2To] that extend the 
previous one. We iterate this construction until covering the interval [0,T]. □ 

3.6 Passing to the limit e — )■ 

We construct a solution of the non-regularized Nanowire drift-diffusion Poisson system by 
passing to the limit e — )■ in the regularization. First, we recall a statement of an Aubin- 
Lions lemma [21 [12] : 

Lemma 3.11. Take T > 0, g G (l,+oo) and let { fn)n£'N be a bounded sequence of functions 
in L''{[0,T]; H) where H is a Banach space. If {fn)n£N is bounded in L''{[0,T];V) where V 
is compactly embedded in H and dfn/dt is bounded in L'^([0,T]; V) uniformly with respect to 
n G N, then, (/n)neN is relatively compact in L'^{[0,T]; H). 

Proof of Theorem II. 6L We fix T > 0. From Proposition 13.101 there exists and 
solution of the regularized system NDDP^ with the initial data Np^. The bound on 
IliYlll^i^^^-) (|3.20p and the dissipation estimate (Corollarv 13. 6|) furnish a bound of ^^/N^ in 
L°°{{0,T];L'^{uJx)) and in L'^{[0,T]; H^iuj^)). Thus, is bounded uniformly with respect to 
e in L'^{[0,T];W'^'^{u}x)) (since we have the equality d^N^ = l^/j^ldr^y^). Next, using the 
Cauchy-Schwarz inequality and Assumption 11.31 we obtain 

[ ([ \d^Nl + Nld^V^\dxfdt<cf V\t)dt 

Jo ^Juix ' Jo 

where V' is the entropy dissipation rate defined in (|3.18p which is bounded in L^([0, T]) 
uniformly with respect to e. From the drift-diffusion equation ()3.7p . we conclude that dtNg 
is bounded in L'^{[0,T];W~^'^ (ujx)) uniformly with respect to e. Therefore, we can apply the 
Aubin Lemma 13.111 for q = 2, H = L^{ujx) and V = W^'^{ujx)- There exists a subsequence 
(that we still denote abusively iV|) such that — ^ Ns strongly in L^([0, T]; L^{ujx)). Finally, 
for this function Ng, we solve the unregularized Nanowire Poisson system and construct V 
such that V E L'^{[0,T]; H^{n)) (Proposition [32]) and lim.^o IIV'' - V\\L^iio,T];m[n)) = 
(thanks to Proposition 13. 9p . 

The last step is to pass to the limit e — )• in the drift-diffusion equation. We have 
rT 

NIdxVgdxdt < C||A''|||ii([o,T];L2{i.j^))ll^^llL°°([0,T];Hi(n))- 
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Corollary 13.61 shows that ||-^|||Li([o,r];L2(ij^)) is bounded independently of e and we conclude 
that there exists a nonnegative constant C independent of e such that 

f f Nld^V^dxdt < C. (3.30) 

It gives a sense to the drift-diffusion equation when e — )• 0. Finally, using (|3.4p . we immediately 
deduce that V^^ Vnn in 1?{^,T\]H^{ijj^)) and that 

converges in L^([0,T] x w^:)- It is enough to prove that 

Nid.,yi - iv,a,y, in p'([o, r] x w,). 

Thus, up to an extraction, (Ng^ ^) is a solution of the NDDP system. Moreover, by semicon- 
tinuity, we can pass in the limit in the a priori estimates such that we still have the relative 
entropy estimation of Proposition [331 for {Ng^V). □ 

Appendix : Technical Lemmas for the diffusive hmit 

Lemma A.l. For all function p i— t- ^{p) polynomially increasing as well as all its derivative 
and for all (p G C°°([0, T] x M) compactly supported, the function if) = 70 can he taken as test 
function in the weak formulation \2.2'T\ j of the Boltzmann equation. 

Proof. Let p ^ (,r{p) such that G V{[-R,R]), < < 1, |5p^ij| < 1 and 1 a.e. 

when R — t- +oo. We set tpR = (pj^B.^ function with which we can write the weak formulation 
()2.27p . To pass to the limit R — +00, it suffices from a Lebesgue theorem that 7/n G Lp{R) 
and p^fn G Lp{^) as well as for dp^. However, with the Cauchy-Schwarz inequality, 

+ \p\)dp <( f{l + \p\?Mr.{p)j\p)dp] ( [ ^{^dp] < 00, 



M„(p) 

because 7 is polynomially increasing. □ 

Lemma A. 2. Let Q be defined in i2.8\} . There exist nonnegative constants Cq, Ci and C2 
such that V(t, x) G [0, T] x M .- 



Ci{l + \p\)< 



<C2(1 + H), (A.l) 



Mn 

< Co(l + IpI^) for s=t, X and p. (A.2) 



M 

Proof. By the definition ()2.8p . we have 

KQn = Qt{&)n + -^Mn- (A.3) 

where we denote Q'j^{Q)n = -A^nl]n'>i/ (^n,n'Qn'dp' and A„ = I]n'>i / ctn,n'Mn'dp' . Using 
Assumption 12. H we immediately find ai < < a2- Applying a Cauchy-Schwarz inequality. 



|Qj(e)„| < a,Mn 2 / \Qn'{p')\dp' < a2Mn(^ I ^^^fl dp' 

n'=l \ n'=l / 

Since G Lj^, (jA.ip follows directly from ()A.3p and we differentiate it to obtain ()A.2p . □ 



2 \ 1/2 
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